Lane formation in driven colloidal mixtures: is it continuous or discontinuous? 
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Binary mixtures of oppositely charged colloids driven by an electric field are shown to exhibit 
a nonequilibrium transition towards lane formation if the driving force is increased. Brownian 
dynamics computer simulations and real-space experiments are employed to study hysteresis effects 
in an order parameter measuring the extent of lane formation upon increasing and decreasing the 
driving force. Both from simulation and experiment, we find that lane formation due to electrical 
fields is continuous. However, simulations show a discontinuous transition if the driving force is 
gravity. 
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In equilibrium, there is a fundamental difference be- 
tween a sharp phase transition which exhibits a jump 
in a certain derivative of the free energy with respect 
to a thermodynamic variable [1] and a continuous cross- 
over where no such discontinuity exists in the thermody- 
namic limit. This is much less clear for nonequilibrium 
phase transitions since a free energy does not exist in 
general in nonequilibrium. Here, in many situations, an 
order parameter can still be defined in the steady state 
and hysteresis behaviour typically serves as a criterion to 
discriminate between discontinuous and continuous be- 
haviour Q. While it is by now well- understood how the 
order of an equilibrium phase transition is controllable 
by the interparticle interactions [H, Q], the question of 
which key parameters determine the existence and order 
of nonequilibrium phase transitions is much less under- 
stood due to the multitude of parameters character- 
izing the nature of the dynamics which do not apply in 
equilibrium. 

Lane formation in a binary mixture of Brownian par- 
ticles which are driven by a constant external force de- 
pending on the particle species [f| represents a prototype 
of a nonequilibrium transition in a continuous (i.e. off- 
lattice) system [f|. Brownian dynamics computer simu- 
lations strongly support the scenario that for increasing 
driving force the system undergoes a transition from a 
mixed steady state towards a steady state where macro- 
scopic lanes are formed. Such a laning transition has 
been found in models for oppositely driven repulsive mix- 
tures [1, 0, H, M, EH EH in two and three spatial dimen- 
sions. In two spatial dimensions, it was found [f| that 
a suitable order parameter which detects laning exhibits 
a significant hysteresis - if the driving force is increased 
and subsequently decreased - which signals a discontin- 
uous nonequilibrium phase transition. This issue is still 
unclear in three spatial dimensions. Though thegen- 
eral scenario occurs for pedestrian dynamics [H, Il3l ] in 
driven granular 3, 15| and colloidal [n| matter and in 



complex plasmas |17| , a quantitative comparison of lan- 
ing with simulation data and the determination of the 
underlying order of the transition is still lacking. 

The aim of this letter is twofold: first, we show that op- 
positely charged colloidal particles driven in an electric 
field [161 ] exhibit lane formation in quantitative agree- 
ment with our Brownian dynamics computer simula- 
tions. Second, we address the question whether the 
transition towards lane formation is discontinuous or a 
smooth crossover. While charged colloids are found to ex- 
hibit lanes in a continuous way upon increasing the elec- 
tric field strength, additional simulations reveal that the 
transition becomes discontinuous if gravity is the driv- 
ing force. The physical reason is that hydrodynamic in- 
teractions are screened in the electric context [3, EH 
while they are long-ranged for gravity. In fact, three- 
dimensional simulations which neglect hydrodynamic in- 
teractions completely yield a continuous transition as 
well. Our results demonstrate that the existence and or- 
der of a nonequilibrium phase transition depends on de- 
tails of the dynamics even when the particle interactions 
which entirely determine the order of phase transitions 
in equilibrium are kept fixed. 

In our experiments, we used suspensions of poly- 
methylmethacrylate colloids in a density matching mix- 
ture of cyclohexyl bromide and cis decalin into which 
we dissolved 60 /xM tetra butyl ammonium bromide 
salt. An equimolar binary suspension of colloids fluo- 
rescently labelled with Rhodamine and 7-nitrobenzo-2- 
oxa-l,3-diazol with volume fraction t)n — "Or = 0.1 was 
prepared. The diameter of both species was 1.2 fim. The 
inclusion of salt led to both colloidal species accquiring 
opposite charges of Zi = lOOe, and a Debye screening 
length of k^ 1 w 300 nm |l6[. The solvent dielectric con- 
stant e was 5.5, and shear viscosity was rj = 2.2 cP. 
All experiments were conducted at room temperature 
T = 300 K. We used a Leica NT confocal laser scanning 
microscope. Experimental data was taken in 2D scans of 
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typically 100 fim in width at a frame rate of around 2 
frames per second, ie coordinates were collected from 
a diffraction-limited slice of ~ 2/ito in depth. 

In our Brownian dynamics computer simulations, we 
consider an equimolar binary mixture of N = 1024 op- 
positely charged colloidal particles of hard-core diameter 
a at corresponding total volume fraction <fi = 0.2 dis- 
persed in a solvent whose viscosity and dielectric con- 
stant matched the experimental values, likewise we set 
and T = 300 K. 

The simulation is performed in a cubic box of length I 
with periodic boundary conditions in all three directions 
[20| . Apart from their steric repulsion, the particles in- 
teract via the screened Coulomb potential 



Vij(r) 



,(-K<r(r/cr-l)) 



(1) 



e(l + Kcr/2) 2 r 

where i and j label the particle species, Zi denotes the 
particle charge, r is the interparticle separation and k is 
the inverse Debye-Huckel screening length depending on 
the salt concentration in the solvent. We take k<t = 4 
in order to match the experimental salt concentration. 
An electric field E yields the driving forces F° xt = Z*E 
acting on the particles. The actual electrophoretic charge 
renormalization is taken to be Z* /Zi = 0.8 as obtained 
from the averaged drift velocity in the electric field at 
low densities [2 lj . 

Time-dependent trajectories of the particles are 
calculated using a finite-time step method with a 
configuration-dependent diffusion tensor 22j. The mag- 
nitude of friction is set by the Stokes-drag expres- 
sion where the free diffusion constant is given by 
D^ksT I "iixr\a h (with ks denoting Boltzmann's con- 
stant) which sets the typical Brownian time scale tb = 
Do/a 2 . Here the hydrodynamic diameter of the parti- 
cles is set very close to the interaction core a such that 
<jh = 0.98(7. We use the pairwise approximation based 
the Long-Ajdari mobility tensor in order to include 
hydrodynamic interactions properly for a driving electric 
field [19j. After an initial relaxation period of typically 
20tb the system runs into a steady state. The time-step 
used in the simulation was At = 10~ 4 T£. 

Experimental snapshots in the steady state for two 
electric field strengths E — 30kV/m and E — lOOkV/m 
are presented in Figure 1. The vertical extent of the op- 
tical slice is set by diffraction, leading to a lengthscale 
of order 2/jto. The snapshots clearly reveal the onset of 
laning for increasing electric driving field. Corresponding 
snapshots of our Brownian dynamics computer simula- 
tions are presented in Figure 2. For sake of comparison, 
the same field strengths were used as for the experimental 
snapshots. The simulation configurations show a similar 
tendency towards laning along the drive direction when 
the electric field is increased. 

In order to quantify the extent of laning, we define a 
laning order parameter m as follows: we assign a cylin- 





FIG. 1: Typical particle configuration for positively (green 
spheres) and negatively (red spheres) charged colloids in an 
electric field of strength E along the horizontal direction. 
Only particles whose centres are within a slice of thickness 
2jj,m are shown. The length bar is 10 /im. The formation 
of lanes can clearly be seen, (a) experimental snapshot for 
E = 30kV/m, (b) experimental snapshot for E = WOkV/m. 
The other parameters are given in the text. 



der of diameter 0.9a and total height 6<r symmetrically 
around the center of each particle i such that the cylin- 
der axis is parallel to the driving field. If one or more 
particles of different species are contained in this cylin- 
ders, we define the value of m, = to this particle while 
777-2 — 1 if the cylinder is free from other particles with 
an opposite charge. The averaged dimensionless laning 
order parameter 



N 

m =< rrii 
i=i 



> /N 



(2) 



where < ... > refers to a particle and steady-state av- 
erage and N is the total number of particles considered 
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FIG. 2: Computer simulation snapshots for positively (green 
spheres) and negatively (red spheres) charged colloids in an 
electric field of strength E along the horizontal direction. 
Only particles whose centres are within a slice of thickness 
2a are shown, (a) E = 30fcV/m, (b) E = WOkV/m. The 
other parameters are given in the text. 



measures the averaged extent of laning. By construc- 
tion, the laning order parameter m equals 1 if particles 
of the same species are always on top of each other along 
the drive but it vanishes in a completely mixed situation. 
We have averaged this laning order parameter over many 
simulation configurations. For the experimental data, a 
similar laning order parameter was obtained by using a 
projected cylinder (i.e. a rectangular cut) to the plane 
containing the particles. Plots of the laning order param- 
eter m versus the strengths of the electric field applied 
are presented in Figure 3. 

First, there is reasonable agreement between simula- 
tion and experiment. The increase of the laning order 
parameter with the external drive is smooth. The ex- 
perimental data reveal a slightly smaller laning order pa- 
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FIG. 3: Averaged laning order parameter m versus elec- 
tric field strength E (in units of kV/m). The open squares 
(full line) are the simulation results upon increasing the field 
strength, the full triangles (dashed line) shows the same upon 
decreasing the field strength. No hysteresis is found. The 
diamonds (dot-dashed line) are the experimental data. The 
parameters are the same as in Fig. 1. 



ramcter than the simulation which is acceptable given 
the uncertaincy in the particle charge. The experimen- 
tal data were taken both increasing and decreasing the 
external field, and no difference was observed. Hence we 
conclude that effects of hysteresis are largely absent. We 
remark that possible effects arising from electro-osmotic 
flow are not exactly known. 

Finally, we consider the case of gravity in which the 
driving forces result from a different buoyant mass of the 
two particle species. In particular, we treat the symmet- 
ric case here where plus and minus charged particles feel 
the same gravitational force but just their sign is differ- 
ent. Our Brownian dynamics computer simulations are 
now performed with the unscreened Rotne-Prager mo- 
bility tensor [23| replacing the screened Long-Ajdari ten- 
sor. All the other parameters are kept unchanged. In 
Figure 4 the simulation results for the laning order pa- 
rameter are shown for increasing and decreasing driving 
forces. The results are compared to the case of an electric 
field where the mobility tensor is screened due to countc- 
rion counter-motion. A significant hysteresis behaviour 
is detected for sedimentation which clearly signals a dis- 
continuous transition towards laning while no such hys- 
teresis is present for electrophoresis and in the complete 
absence of hydrodynamic interactions. This clearly indi- 
cates that it is the range of hydrodynamic interactions 
which make lane formation discontinuous. In a strictly 
two dimensional system, laning occurs in a discontinuous 
way even when hydrodynamic interactions are ignored 
completely Q . Therefore besides the hydrodynamics the 
system dimensionality plays an important role in deter- 
mining whether laning is discontinuous or continuous. 
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FIG. 4: Averaged laning order parameter versus reduced 
driving force /* both upon increasing (open squares, full 
line) and decreasing (full triangles, dashed line) the driving 
forces: (a) for an electric driving field where /* = E mod Z* 
mod a jknT and hydrodynamic interactions are screened, (b) 
for a gravitational force Mg where /* = ±Mgo /ksT and hy- 
drodynamic interactions are long-ranged, (c) same as (a)/(b) 
but for neglected hydrodynamic interactions. The length of 
the cylinder in the definition of the order parameter here spans 
the whole simulation box. 



In conclusion, we have shown that oppositely charged 
colloids driven by an external electric field show a con- 
tinuous tendency towards lane formation if the driving 



strength is increased. Real-space experimental data are 
in quantitative agreement with Brownian dynamics com- 
puter simulations which include the screened hydrody- 
namic interactions between the driven colloids. The con- 
tinuous crossover towards laning is in marked contrast 
to sedimenting colloids where lane formation occurs via 
a nonequilibrium first-order phase transition as revealed 
by a significant hysteresis behaviour in a suitable order 
parameter. This change must be attributed to the long- 
range nature of hydrodynamic interactions in the sedi- 
mentation case. In a reduced system dimensionality, on 
the other hand, laning is discontinuous as well even when 
hydrodynamic interactions are ignored. Hence the order 
of the lane formation is determined by both, hydrody- 
namic interactions and system dimensionality. Therefore 
the more general conclusion is that the existence and or- 
der of a nonequilibrium phase transition depends on de- 
tails of the dynamics even when the particle interactions 
which entirely determine the order of phase transitions 
in equilibrium are kept fixed. 

Lane formation does not only occur in colloidal systems 
but also in dusty plasmas where hydrodynamic interac- 
tions are absent but inertia effects play a dominant role 
[l7i | . We expect that, in three dimensions, the tendency 
towards lane formation will be continuous there as well. 
In strict two dimensions, e.g. for pedestrian dynamics 
and shaken granular matter, we expect a continuous be- 
haviour irrespective to the details of the dynamics. 
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